Measuring the Entanglement between Double Quantum Dot Charge Qubits 
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We present a scheme for creating and measuring entanglement between two double quantum dot 
charge qubits in a transport set-up in which voltage pulses can modify system parameters. Detec- 
tion of entanglement is performed via the construction of a Bell inequality with current correlation 
measurements. An essential feature is the use of the internal dynamics of the qubits as the con- 
stituent electrons tunnel into the leads to give the single-particle rotations necessary for the Bell 
measurement. 

PACS numbers: 03.65.Ud, 73.63.Kv, 73.50.Td, 73.23. Hk 



In an important recent experiment Shinkai et al. 
have demonstrated correlated coherent oscillations be- 
tween two coupled double quantum dot (DQD) charge 
qubits formed in a top-gated semiconductor heterostruc- 
ture 01 ■ From transport measurements, indications were 
obtained that it should be possible to perform a suite 
of universal two-qubit quantum gates with such a setup. 
However, such a claim can only be substantiated if it can 
be shown that the operations can entangle the qubits Q . 
This brings us to the question that is the focus of this 
Letter: If the qubits of Shinkai et al were entangled, how 
could we tell? Is it possible to detect and measure the 
entanglement between DQD charge qubits in a transport 
set-up such as that of Ref. [ij? 

We answer these questions here by describing a series of 
shotnoise measurements that can be used to construct a 
Bell's inequality (BI) Q , the violation of which provides a 
clear signal of, and in certain circumstances quantitative 
information about, the entanglement between the qubits. 
Shotnoise and the Bell's inequality have been combined 
to study entanglement in mesoscopic systems before Q. 
However, the system here differs in several respects. Most 
important is that here we are in the sequential-tunnelling 
regime and coherence between the electrons is assumed 
to be lost once the electrons tunnel to the leads. Further- 
more, the BI requires single-qubit rotations and it is not 
initially obvious how this may be accomplished. In this 
Letter, we show how these obstacles can be overcome 
by making use of the internal dynamics of the qubits. 
As the qubit electrons tunnel into the leads, i.e. as the 
qubits decay, they experience the action of the system 
Hamiltonian and this rotates the qubits. We show how 
an appropriate set of current correlation measurements 
can extract the relevant information from the stochastic 
background of the qubit decay. 

The set-up of interest is sketched in Fig. [T^ with four 
quantum dots (QDs) making up two charge qubits, with 
the position of the electron within a DQD (left or right) 
corresponding to the two logical qubit states. With a 
pseudo-spin convention that | j) = \L) and | f) = \R), 
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FIG. 1: (a) The system consists of two DQDs with each 
dot connected to its respective lead. With one electron in 
each DQD, two charge qubits are formed that interact with 
interaction strength J. (b) The maximum concurrence C of 
the pure two-qubit state evolving from state \LL) a.t t — 
under the action of Hoqu- Only the time interval [0, Tmax] 
with Tmax = IS Considered. Other parameters were e^^-* = 
= e and T^^) = T^-^) = T. 



the Hamiltonian of the isolated two-qubit system is 
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with ei and Ti the detuning and tunnel coupling of DQD 
i = 1,2, and with J the strength of the Coulomb inter- 
action between the electrons. 

We assume that the DQD parameters as well as the 
chemical potentials of the leads can be controlled to a 
limited extent via voltage pulses and, as in Refs. [J, 0], 
we envisage an experiment that consists of a series of re- 
peated steps. At the start of the cycle, both DQDs are 
empty. Then the left chemical potentials are raised, such 
that one electron tunnels into each of the dots from the 
left. As in Refs. maintaining a large bias across the 

dots produces detunings Ci such that the electrons tunnel 
into localised states li*-*-*). The right chemical potentials 
are then raised such that the electrons are trapped and 
two qubits are formed. They remain trapped in the sys- 
tem with the full two-qubit Hamiltonian acting on them 
for a time Tjnit , at the end of which the two-qubit system 
is left in the state po. We want to measure the entan- 
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glement of this state. Finally, all chemical potentials are 
dropped below the dot levels and the electrons escape 
into the leads. Sufficient time is waited for this to have 
happened with certainty before the sequence is repeated. 

The only control over dot parameters assumed here is 
that they be set once during the initialisation of state 
po, and subsequently to a different set for the detection 
phase. It is essential that various different configurations 
are possible so that there is a range of single-particle rota- 
tions for the BI measurement. As we show below, with an 
appropriate choice of current correlation measurements, 
we can construct the correlator f (a, b) = (a-cr^^^b-cr'^') 
with unit vectors a and b determined by the DQD pa- 
rameters. Repeating the experiment with different pa- 
rameters allows us to build the correlator 



F = £{a, b) + £{a, b') + £{a', b) - £{a' , b'). 



(2) 



The pertinent BI is the CHSH inequality, which reads 
-2<i^<2[Z]. A measurement of |F| > 2 therefore 
indicates entanglement. Maximising \F\ with respect to 
vectors a,b,a',b' yields |Fmax|, and for pure states we 
have |Fmax| = 2vTTC^ with C the concurrence, a 
measure of two-qubit entanglement For mixed states, 
the region of violation is bounded by max [2,2^/2C] < 



C2 



li^maxl < 2VT 

Initialisation. We first discuss the nature of the en- 
tangled states that can be produced within this scheme 
before analysing the detection protocol in more detail. 
The Hamiltonian of Eq. ([T]) can be used to form many 
different quantum gates. For example, with ei = —£2, 
Ti ~ T2, and |J ± ei| ^ Ti, the Hamiltonian is 
effectively that of a FLIP gate, ffDQD ~ ^^flip = 
if7 {\LL){RR\ + mR){LL\), with effective coupling O = 
2Tf J/(j2 — el) Starting with the system in state 
|LL), and allowing it to evolve for a time t, the density 
matrix of the system becomes 



p(t) = lj{l + a)\LL){LL\ + {l-a)\RR){RR\ 
+iCi\LL){RR\ - \RR){LL\)} , 



(3) 



with a ~ cos(r2r) and C = sin{QT)e~'^^ , the concurrence 
of the state, in which we have included a dephasing at 
rate 7. These states, although mixed for finite 7, have 
the same violation of the Bell inequality as a pure state 
of the same concurrence, F{p) = 2\/l + [ll|. The 
entangled states obtained by letting Huqd act on the 
state \LL) for a time arc not restricted to these FLIP- 
class states. The entanglement of this more general class 
of states is investigated in Fig. [TJ^. Clearly, a high de- 
gree of entanglement (concurrence approaching unity) is 
obtainable under reasonable conditions. 

Entanglement detection. Let us initially consider an 
ideal model of the decay of the two-qubit system. We 
assume that, once the initial entangled state is created, 
the interaction between the qubits is turned off (J = 0). 



This will be a good approximation if the coupling of the 
QDs to the leads is strong enough that the 'measurement 
phase' is quick compared to the interaction timescale 
J~^. Let us also assume the most general form for the 
single-qubit Hamiltonian H^'^'> = A'^'^n^*' -cr^'), where A*^'^ 
is the single-qubit energy scale, n^*) is a unit vector, and 
(t(*) is a vector of Pauli matrices. We return to our more 
realistic DQD model shortly. 

We consider the DQD system to be in the strong 
Coulomb blockade regime, such that at most one excess 
electron is permitted in each DQD. We describe the tun- 
nelling of electrons with a Markovian master equation, 
p = Cp, with p the density matrix of the qubit pair and C 
the total system Liouvillian. The initial state of the sys- 
tem is po, some entangled state such as that of Eq. Q. In 
the noninteracting approximation (J = 0), we can anal- 
yse the dynamics of each qubit separately and need only 
consider the two qubits together when we introduce the 
initial state. To facilitate calculation of current statis- 
tics, we introduce counting fields Xl^ and of elec- 
trons in the L and R leads of dot z = 1 , 2 to obtain 
the x-resolved master equation: p^'^ (x) = (x)p^*^ ix) 
with £('^(x) = ^^^Hxl^ 1X11)1 x-resolved LiouviUian 
of dot i, and similarly for the density matrix, which in- 
cludes the empty state jO^*)) as well as the two qubit 
states |L(*)) and The Liouvillian can be written 



£W(X) 



E('^(x) with first term describing the 



internal dynamics Lq ^ p^'^^ ~ —i [i/^^^p^*)], and the sec- 
ond, the coupling to the leads. With chemical potentials 
set far below the dot levels, the latter can be written 



S]W(X)P 



(0 



4 E {^f^^p^ 



a=L,R 



},(4) 



with operators = |0)(a^'^| describing the jump on an 
electron from localised state a^*^ = L'*^ i?'*^ to the leads. 
Setting Xa = in Eq. ^ we obtain the familiar master 
equation in Lindblad form which describes the evolution 
of the actual density matrix p*-*'. From this point on, 



we assume that all rates are identical, F 



(i) 



r 



r, 



and assume that this rate is faster than dephasing rate 
7, such that such dephasing from external sources can be 
neglected in the detection phase (see later). 

In Laplace space, the density matrix of the system 
at arbitrary time is p*^*^(x;z) = n^'^\xi^)Po with the 
X-resolved propagator for a single DQD, ri(*)(x;z) = 



sW(x) 



We are only interested in the 
situation in which the system starts with one electron 
in each DQD, in which case, the propagator in the long 
time limit reduces to 



J, 



+ J, 



— +7^(0«,n«)) .(5) 



^a'^ are jump super-operators defined 
by ji^V^'^ = |0«)(a«|p«|a«)(0«| and 
7^(6'W,n(*))pW = [/(*)pWj7W^ with U'^^ 
exp (— |(?'*-'n'^*' • cr^*') a unitary rotation about axis 
n(') by angle 6l(*) = \ arctan (2A(*) /F) . The above form 
of the propagator is the main formal result of this work. 
It shows that, in the long-time limit, the behaviour of 
the system effectively decomposes into two parts: one in 
which the qubits decay directly, and one in which the 
qubits are first rotated and then leave the dots. This 
rotation originates from the action of the single-qubit 
Hamiltonian i/*-*' acting for a time governed by the ratio 
of A(*) to r. 

The moment generating function (MGF) for the 
two-qubit system in the long-time limit is M{x) ~ 

Tr |riTC''(x) (g) r2^^(x)po|, where po is the two-qubit 

entangled state. The first moment, M\ = 

dM I d{ixx) , corresponds to the mean number of 

electrons transferred to lead X = L,R oi dot i in the 
measurement part of the cycle. More importantly, let us 

as the cross 



define MxY = ^^^(/^(ixx )^(«Xy^) 



correlator between the number of electrons emitted into 
lead X of dot 1 and into lead Y of dot 2. This quantity 
can be extracted from shotnoise measurements [5| as the 
preparation-detection cycle is repeated continuously. In 
analogy to the standard CHSH measurement 0], let us 
define 



C(7^(l^7^(2)) = Mll-Mlb-Mrl+Mrr 

X (l(2)+7^(0(2),n(2)))pQ|. (6) 



This quantity still has contributions from the non- 
rotating decay. However, if we take the following combi- 
nation of C-correlators: 

£(a,6) = 4C(7^(l^7^(2)) -2C(7^(l^l(2)) 

-2C(l(l^7^(2)) + C(l(l^l(2)), (7) 

we obtain £{a,b) = {a ■ cr'^-^^b ■ cr'-^^) with a ■ er^^' = 

and similarly for b. This series of mea- 
surements then yields exactly the correlation function 
required to form the CHSH inequality. The operations 
can be realised simply by any rotation about 
the z-axis, i.e. with ^ Ti. 

Since we start with exactly one electron in each DQD, 
in the long time limit the total number of electrons to 
leave each DQD must also be unity. This means that 
we can rc-express the correlation function of Eq. ^ as 
C(7^^l^7^(2)) = 1 + 4Mbr - 2(X^^' This has 

the great advantage that one need only measure currents 
on one side of the dots. Furthermore, the electron count- 




FIG. 2: (a) Time evolution of CHSH correlator F for two 
initial FLIP-class states (Eq. (|3}): one separable (C — 0, 
dashed lines) and one maximally entangled {C — 1, contin- 
uous lines) with J — 0. The thick lines show F obtained with 
vectors close to those yielding maximum CHSH violation for 
C = 1. The thin lines show F with a choice of vectors such 
that the inequality is never violated. The violation (or other- 
wise) becomes apparent after a time ~ 4r~^. (b),(c) Plot of 
|Jincas| vs |Fniaxi for FLIP-class statcs, with |-Fmcas| obtained 
the under constraint < 9 < Smax, and with finite interaction 
strength J. The straight line indicates |-Fmeas| = |i^max| 



ing we have pursued here only accounts for electrons leav- 
ing the dots and ignores those entering the dots in the ini- 
tialisation phase. By expressing the measurement purely 
in terms of right-lead quantities, we avoid having to ex- 
plicitly take the latter into account. Figure [5^ shows the 
CHSH correlator, F, as a function of time. Violations of 
the CHSH inequality become visible after a time ~ 4F~^ 
Non-idealities. The situation in real dots differs from 
the preceding analysis in two respects that we now ad- 
dress: (i) the single-qubit Hamiltonians do not have a ay 
component and are further restricted to experimentally 
accessible DQD parameters; and (ii) it is unlikely that 
the interaction between qubits can be completely sup- 
pressed during the measurement phase. The single-qubit 
part of -ffoQD can be rewritten as 



A(') (sin(i0W)(Ti''+cos(i0('))ai') 



(8) 



with A(^) = Vej+T^- If and e, are unrestricted in 
magnitude, angle (/)*^*) has the range — tt < 0^'-' < tt, and 
angle S*^*' is bounded as < 6*^'^ < tt with the upper 
bound attainable only in the A^^'/F oo limit. With 
these bounds, a sufficient range of rotations can be per- 
formed to ensure maximum violation, despite the absence 
of the (7 J, component. However, since A('Vr is finite, 
the rotation angle is restricted, < 9^''^ < 6'max, with 
^'max = ^ arctan (2A(*) /F) . This constraint limits the 
optimisation of F. More importantly, with interaction 
during the measurement, the value of F obtained with 
the above scheme is no longer equal to the true value 
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^max = 0.997r 




FIG. 3: As Fig.[2]3, but for randomly chosen density matrices. 
Ideal values of I-Fmaxl calculated from Refs. 



density i 



of F of the state. Figures [2] and [3] show the maximum 
value of \F\ that would be measured, |Fmcas|, against 
the value of |Fmax| calculated from the state. The mea- 
sured results arc obtained from numerical integration of 
the Master equation with the interacting Hamiltonian of 
Eq. H]) and subsequent optimisation of F subject to the 
above constraint on 6^'^'' . 

Figures [Ud and c shows the results for initial states 
belonging to the FLIP class. In the noninteracting case 
(J = 0) the restriction to 6'max < ti" causes only very small 
reductions in the measured value of |-Fmax| provided that 
e^i^/iT > 0.7, which corresponds to (A(*)/r)^^^^^ > 1. 
Furthermore, it is clear from the finite- J results that, for 



FLIP-class states, \F„ 



< \Fr 



even with interaction 



present, and thus |i^mcas|, and the concurrence calculated 
from it, provide lower bounds for the actual values of 
l^inaxl and C — the measurement never flags an unentan- 
gled FLIP-class state as an entangled one. Note further 
that the larger ^max is, the smaller the deviations. For ex- 
ample, at ^max/"' = 0.99, even a value of J = F produces 
only a deviation in iFmaxI of ~ 5%. Figure [3] shows the 
same comparison for a set of randomly generated density 
matrices. In this case, the measured value I-Fmoasl is no 
longer a lower bound for |i^max|- For small J, however, 
the measured values are strongly clustered around the 
ideal values, and tight upper and lower bounds can be 
derived. For larger values of J, the clustering is not as 
tight, and unentangled states can give quite large values 
of l^lnoasl- From these data, we infer that a reliable signa- 
ture of entanglement is |-Fmoas| ^ 2.35 for ^max = 0.997r 
and |F,ncas| ^ 2.2 for 6'max = O.Ttt. From this point of 
view, tighter restriction of ^max niay be advantageous. 

From the foregoing, we can identify the parameter re- 
quirements for obtaining the strongest possible signal of 
entanglement. In order to obtain tight bounds on -Fmax, 



we require J < F < Amix in the measurement phase. 
Obtaining a high entanglement in the initialisation phase 
implies a different relation between A*^*) and J — from 
Fig. [2] we observe that high concurrence occurs when 
Ti ~ J. If we assume Amlx ~ (T^Omax, this implies that 
the single-qubit coupling should change from Ti ^ J to 
Ti > J between phases. This may be feasible, but note 
that this result was obtained with a maximum initiali- 
sation time of J~^. If we allow an initialisation time of 
10 J"^, say, then due to the greater range of oscillation 
explored, high concurrences require Ti ^ A^'^ ^ O.IJ, 
such that Ti > J in both phases. In the end, the maxi- 
mum initialisation time is determined by the dephasing 
rate 7 and we require 7 J for successful operation. For 
a dephasing time of Ins Q (corresponding to 7 ~ l//eV), 
and interaction strength J = 25^eV, this relationship is 
satisfied. This also implies that F 2> 7, consistent with 
the neglect of debasing during the read-out phase. Note 
that it might also be experimentally possible to increase 
the potential barrier between the two DQDs at the end 
of the initialisation phase, thus separating the electrons 
and reducing their interaction. This would further im- 
prove operation. 

In summary, we have described a way of obtain- 
ing information on the entanglement of DQD charge 
qubits which uses shotnoise measurements of the decay- 
ing qubits to determine the CHSH parameter i^max- This 
BI approach does not need complete control over single- 
qubit rotations and it is not even necessary to know which 
rotations were performed. This can be contrasted with 
density matrix tomography IJ], which requires precise 
control and knowledge of the single-qubit rotations. Our 
approach might also be used to measure DQD entangle- 
ment generated in other ways, e.g. through interaction 
with a common bath Finally, we mention that 

the use of internal dynamics to generate BI rotations may 
have broader applicability to the study of entanglement 
in other decaying system, e.g. quantum optics. 
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